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Abstract 

Our goal is the mathematical analysis of a two phase (liquid and gas) two 
components (water and hydrogen) system modeling the hydrogen displace- 
ment in a storage site for radioactive waste. We suppose that the water 
is only in the liquid phase and is incompressible. The hydrogen in the gas 
phase is supposed compressible and could be dissolved into the water with 
the Henry's law. The flow is described by the conservation of the mass of each 
components. The model is treated without simplified assumptions on the gas 
density. This model is degenerated due to vanishing terms. We establish an 
existence result for the nonlinear degenerate parabolic system based on new 
energy estimate on pressures. 

Keywords: Degenerate system, nonlinear parabolic system, compressible 
flow, porous media0. 



1. Introduction 

An important quantity of hydrogen can be produced by corrosion of fer- 
rous materials in a storage site for radioactive waste. A direct consequence 
of this production is the growth of hydrogen pressure around alveolus. This 
increasing gas pressure could break the surrounding host rock and fractures 
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could appear in the confinement materials. This problem renews the mathe- 
matical interest in the equation describing multiphase/multicomponent flows 
through porous media. The cases of immiscible and incompressible flows have 
been treated with "global pressure" introduced by G. Chavent and J. Jaffre 
[sj by many authors, we refer for instance to [HI, 0, 0, [l2| where existence 
results are obtained under various assumptions on physical data. For two 
immiscible compressible flows without exchange between the phase, we refer 
to 14, 15 where the authors obtain the existence of solution when the densi- 



ties depend on the global pressure and to [l7|, [l8| for the general case where 
the density of each phase depends on its own pressure. This approach is also 
used in [2|, |3| to treat a homogenization problem of immiscible compressible 
water-gas flow in porous media. For miscible and compressible flow, we refer 
to 0, llOl for more details. 

In j^, the authors derive a compositional model of compressible multi- 
phase flow in porous media. They focus their study on models where the 
fluid is a mixture of two components: water (mostly liquid) and hydrogen 
(mostly gas). An existence result has been shown in |20| for this model under 
the assumptions of non degeneracy and of strictly positive saturation. 

Recently in (sf, the authors studied a new model of two compressible and 
partially miscible phase flow in porous media, applied to gas migration in an 
underground nuclear waste repository in the case where the velocity of the 
mass exchange between dissolved hydrogen and hydrogen in the gas phase is 
supposed flnite. 

Let us state the physical model used in this paper. We consider herein 
a porous medium saturated with a fluid composed of two phases (liquid and 
gas) and a mixture of two components (water and hydrogen) studied in (l9| . 



As reported in [19|, the author establish the existence of a weak solution, 
under non degeneracy and slow oscillation assumptions on the diagonal co- 
efficients and with small data for the hydrogen. Our aim is to show global 
solution for the degenerate system without restriction on data. 
The water is supposed only in the liquid phase (no vapor of water due to 
evaporation). In order to define the physical model, we write the mass con- 
servation of each component 

dti^sip^ + $.,pj) + div(pfV, + pJV,) - div(p,Df VXf ) = r,(l.l) 
dt{^Sipr) + dw{prYi)=r^. (1.2) 

Here the subscript / and g represent respectively the liquid phase and the gas 
phase. Quantities <D, pf , pj, Pc. = pl + pZ, s^, = pf/pi {Xj' + = l) 
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and represent respectively the porosity of the medium, the density of 
dissolved hydrogen, the density of the hydrogen in the gas phase, the density 
of the a phase {a = l,g), the saturation of the a phase [si + Sg = 1), the 
mass fraction of the hydrogen in the liquid phase, the diffusion- dispersion 
tensor of the hydrogen in the liquid phase. The velocity of each fluid is 
given by the Darcy's law 

V„ = _K^^iiM (Vp, - pMs) , (1.3) 

where K is the intrinsic permeability tensor of the porous medium, kr^ the 
relative permeabihty of the a phase, /Xq, the constant a-phase's viscosity, pa 
the a-phase's pressure and g the gravity. For detailed presentation of the 



model we refer to the presentation of the benchmark Couplex-Gaz [21| and 
0, 20|. The capillary pressure law is defined as 



PciSl) =Pg-pi, 

is decreasing, (^(s/) < 0, for all si G [0, 1]) and Pc{l) = 0. 

The system (1 l.ip - fl 1.2p is not complete, to closing the system, we use 
the ideal gas law and the Henry's law 

]\/fh 

pj = —Pg, pf = M'H%, (1.4) 

where the quantities M^, H^, R and T represent respectively the molar 
mass of hydrogen, the Henry's constant for hydrogen, the universal constant 
of perfect gases and T the temperature. 

By these formulation, the system f l l.ip - fl 1.21) is closed and we choose the 
liquid and gas pressures as unknowns. From fl 1.4p . the henry's law combined 
to the ideal gas law, to obtain that the density of hydrogen gas is proportional 
to the density of hydrogen dissolved 

= Cipf where C, = j^i= 52.51). (1.5) 

Remark that the density of water pf in the liquid phase is constant and from 
the Henry's law, we can write 
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where C2 is a constant equal to H^M^. 

Then the system (I l-ip - d l-2p can be writen as 



dt ($m(sOpf) + div (pfV, + CipfV,) - div (CsX-Df Vp,) = r„(1.6) 



^"^^^ + div (V,) = ^. (1.7) 



where m{si) = si + CiSg. 

Note that the mass exchange between dissolved hydrogen and hydrogen 
in the gas phase is static using the Henry's law opposite then supposed in 
i. 



2. Assumptions and main result 

The main point is to handle a priori estimates on the approximate solu- 
tion. Due to the degeneracy for dissipative terms d\v{p^MaVpa), we can't 
control the discrete gradient of pressure since the mobility of each phase van- 
ishes in the region where the phase is missing. So, we are going to use the 
feature of global pressure to obtain uniform estimates on the gradient of the 
global pressure and the gradient of a capillary term to treat the degeneracy 
of the dissipative terms. Let summarize some useful notations in the sequel. 
We recall the conception of the global pressure as describe in 

M{si)Vp = Mi{si)Vpi + Mg{sg)Vpg, (2.1) 

with the a-phase's mobility and the total mobility are defined by 

M„(s«) = A;,Js„)/p,, M{si) = Mi{si) + Mg{sg). 

This pressure p can be written as 

p = Pg+p{si) =pi+p{si), (2.2) 

with 

dp _ Mi{si) dpc dp _ Mg{sg) dpc 
dsi ~ M{si) dsi dsi ~ M{si) dsi ■ 

We also define the contribution of capillary terms by 

Ml{Sl)Mg{Sg)dp n 

l{s{) = ^^(^^(^0 > and B{s,) = ,{z)dz. 
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We complete the description of the model f l 1.6p - fl 1.7p by introducing bound- 
ary conditions and initial conditions. Let T > be the final time fixed and 
let he n a bounded open subset of M'^ {d > 1). We set Qt = (0,T) x 17, 
= (0, T) X dQ and we note F; the part of the boundary of Q where 
the liquid saturation is imposed to one and r„ = r\r/. The chosen mixed 
boundary conditions on the pressures are 

Pg(t, x) = pi(t, x) = on (0, T) X Ti, 
V, ■ n = V3 ■ n = on (0, T) X r„, 
Xrol'Wpg ■ n = on (0, T) X T,, 

where n is the outward normal to r„. 

The initial conditions are defined on pressures 

p^{t = 0) =pI in Q, for a = l,g. (2.3) 

Next we introduce a classically physically relevant assumptions on the coef- 
ficients of the system. 

(HI) The porosity G W^'^^IQ) and there is two positive constants 0o and 

01 such that 00 < (Pi^) ^ 0i almost everywhere x ^ fl. 
(H2) There exists two positive constants /cq and k^o such that 

||K||(^^(^)),xd < fcoo and (K(x)e,e) > ^olel'-Ve G M^. 

(H3) The functions Mi and Mg G C°([0, 1],M+), M„(s„ = 0) = 0. and there 
is a positive constant mo > such that for all si G [0, 1], 

Mi{si) + Mg{sg) >mo. 

(H4) The densities pa (a=l,g) are in C^(]R), increasing and there exists two 
positive constants pm > and pm > such that 

< pm< PaiPa) < pM, a = l,g. 

(H5) The capillary pressure fonction p^ G C^([0, 1]; M"*") and there exists Pc > 

such that <Pc< |^|. 
(H6) The functions r^^, Vg G L'^{Qt) and r^, Tc, > a.e. for all (t,x) G Qt- 



5 



(H7) The diffusion-dispersion tensor (function of x and si) is a nonlinear 
continuous function of the hquid saturation s; and is bounded for x G O 
and si e [0,1]. In addition, there exist a constant d* > such that 
\fv e M-^, Vx e n, Wsi G [0, 1], (x, si)v, v) > d*\\v\\\ 

(H8) The function 7 G ([0, 1]; M+) satisfies 7(5;) > for < < 1 
and 7(0) = 7(1) = 0. We assume that (the inverse of B{si) = 
'j{z)dz) is an HoldeiH function of order 6^, with < 9 < 1, on [0,i3(l)]. 

Let us define the following Sobolev space 

H^^{n) = {ue H\ny, n = O on Ti} 

this is an Hilbert space with the norm (f^) = || Vn||(i2(Q))d. 

Let us state the main result of this paper 

Theorem 1. Let {H1)-{H8) hold and let the initial conditions (p°,p°) be- 
longs in L^{Q) X L'^{Q) with < < 1. Then there exists a solution {pg,pi) 
satisfying 

Pg G L\0,T;Hl^{n)) and ^WJI^Vpo, G L\Qt). (2.4) 
Si > a.e. m Qt,B{si) G L'^{0,T; H\n)), (2.5) 
^dt{p1{pg)m{si)) G L'{0,T; {Hl^m'), ^Si G L'{0,T; (i/^,(f^))'), (2.6) 
in the sense that for all (f,ip& C^{0, T; (f^)) with (p{T, .) = ip{T, .) = 0, 



^m{si)p1{pg)dt(pdxdt- I <l>m(s?)pf (p°)(^(0, a;)da; 



Qt 



Kpf{pg)Mi{si){Vpi - pi{pi)g) ■ Vifdxdt 
+ Ci / Kp'lipg)Mgisi){Vpg - pg{pg)g) -Vifidxdt 
+ / C2X';"D^Vpg ■ V^dxdt = / rgifdxdt, 

Jqt JQt 

(^sidtipdxdt - / $s°'?/;(0,x)dx 
Jn 

I 'KMi{si){Vpi- pi{pi)g)-Vi^dxdt= [ ^^dxdt, 
Jqt Jqt Pi 



(2.7) 



(2.^ 



^This means that there exists a positive constant h such that for all a, 6 G [0,S(1)], one 
has \B-^{a) - B-^{b)\ < b\a - 
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and the initial conditions are satisfied in the sense that for all ^ G H^^{Q), 
the functions t ^ $p[*(pg)m(s/)^dx, and t J^^SiC,dx, are in C^{[0,T]) . 
Furthermore, we have 

^p1{p,)m{s^)^dx) (0) = / 0pf (p;)m(s°)edx, 
n / Jn 



^si^dx (0) = / ^s^^dx. 
n J Jn 

Remark 1. Remark that, the solutions obtained in the above theorem do not 
satisfy that si < 1, then the functions depend on the saturation are extended 
by continuity for si > 1. 



3. Energy estimates 



The notion of weak solutions is very natural provided that we explain the 
origin of the requirements (1 2.4p - fl 2.6p . In this section, we give estimates on 
the gradient of the global pressure and on the gradient of the capillary term 
B. In order to obtain these estimations, we define gg and Tig by 



dgiPg 



Pa 



Pz'(^) 



dz and Ugipg) := p^{Pg)gg{Pg) - Pg. 



The function Ug verifies n'g{pg) = {p'l^pg))' gg{pg), TigiO) = 0, Tigipg) > 
for all Pg and Tin is sublinear with respect to Pg. This kind of function is 
introduced in 16, 17, 3]. 

By multiplying (1 1.6p by gg{pg) and ( 1 1.7p by Cipi —pg, after integration and 
summation of equations, we deduce the equality 



dt (rn{si)p^{pg)Jgg{pg) + dtSi (Ci pi - Pg 

+ I \iMi{si){\/pi- pi{pi)g)-S/pidx 
Jn 

+ Cl [ KMg{Sg)iVpg - Pg{Pg)g) -VpgdX 

Jn 



dx 



(3.1) 



C2X;'D^Vpg ■ VpgdX 



rgggiPg 



Pt 



-{Cipi-pgUdx. 
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To treat the first term of tlie above equality, let 
M=dt {rn{si)p^{pg)^gg{pg) + dfSi {Ci pi - Pg 



dt im{si)pi{pg)gg{pg)J + i^si{Cipi " j9g ) j + CiSidt{pg - pi) - CidtPg 
dt[m{si)pi {pg)gg{pg)) +dt(si{Cipi -Pg)) +CiSidt{pc) - CidtPg. 



Consider M a primitive of sip'^{si). We can write as = dtS where £ is 
defined by 

£ = m{Sl)Pl{Pg)gg{Pg) + Sl{ClPl - Pg) + ClM{Sl) - ClPg 

= m{Sl) [p^l{Pg)gg{Pg) - Pg) " ClSlPc{Sl) + ClU [Sl] . 



From the definition of the functions Hg and A/", the expression of £ is equiv- 
alent to 

£ = m{si)'Hg{pg) -Ci Pciz)dz. 

Jo 

Integrate f l 3.ip over (0,T), we deduce by using the assumptions {H1)-{H8), 
the positivity of Tig and the sub-linearity of gg{pg), that 

/ Mi\Vpi\^dxdt + / Mg\Vpg\'^dxdt 

Jqt Jqt 

+ [ C2XrD^Vpg-Vpgdxdt<C{l + \\pi\\L^Q,) + \\Pg\\LHQr)), (3.2) 

where C > is constant. In term of global pressure, from the relation fl 2.ip . 
we have the fundamental equality 

M\Vp\' + ^^\Vpf = Mi\Vpi\' + Mg\Vpg\', (3.3) 

The relation f l 2.2p between the global pressure and the pressure of each phase 
prove the following inequality 

WPiWlHQt) + WPgh^QT) < WpWl^Qt) + \\P\\l^(Qt) + Ml^Qt) 



The above inequality and the equality fl 3.3p combined to the estimate fl 3.2p 
ensures that p, Pg G L'^{0,T; H^in)) and B{si) e L'^{0,T; H^{n)). 



4. Construction of a regularized system 

Before establishing Theorem [H we introduce the existence of regularized 
solutions to system fl l.ip - fl 1.2p . First we are interested in a non-degenerate 
system by adding a dissipative term on saturation preserving the positivity 
of the liquid saturation. Precisely, we consider the non-degenerate system: 



(Ci-l)r^div(pfb,)V(p,^-pn) 



r 



9' 



d, - div (KAUs^) Wl - pi{p^>)g)) - V div(V(p^ - p^)) = ^, 

Pi 



completed with the initial conditions (12. 3p , and the following mixed boundary 
conditions, 

p'^g{t,x) =p'^{t,x) = on(0,r)xr, 

(V;' + CiV^ + C^XrD^Vpg + (Ci - l)r/V(p^ - p^)) ■ n = on (0, T) x r„ 

(V;' - r/V(p^ - p"/)) ■ n = on (0, T) X r„ 

(4.1) 

where n is the outward normal to the boundary r„ and = —'KMa{s^)(Vp'2,— 
Paipl)g)- 

We state the existence of solutions of the above system (^r,). 

Theorem 2. Let (H1)-(H8) hold. Letp'^g,p'} e L'^iVt), < s['(x) < 1. Then, 
for all 7] > 0, there exists {Pg,p1) satisfying 

Pl e L\^,T-HlX^)\ $a,(pf(ppm(.n) e L2(0,T; (//^, (!]))'), 
Piip'Mis^) e C° {[0,T]-L\n)),s^ e C {[0,T]-L\n)) , 
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such that for all(p,ip G L'^{0,T; H^^{n)) 



(4.3) 



($a,(m(s;')pf (p^)), ^)+ [ Kpf (p^)M,(sI') (Vpl' - pz(p?)9) ■ V^dxdt 
+ Ci / Kpf (ppM,(.;') (Vp,^ - p,(pP^) ■ Wcpdxdt 

■JQt 

+ / C^iXmD^rVp'l ■ W^dxdt + (Ci - 1)77 / pf (pP V(p^ - p?) ■ V^dxdt 
= / Tgifdxdt, 

Jqt 

(4.2) 

{^dtsl [ KMiis"^) {VpI - pi{p1)g) ■ Vijdxdt 

JQt 

~ V ^{Pg ~ pI) ■ V^/'dxdt = / -^ipdxdt^ 
JQt JQt Pi 

where the bracket (., .) is the duality product between L^(0,T; (ifp^(f2))') and 
L\Q,T-Hl^{n)). 



For the sake of clarity, we omit the index 77 in the problem i^ri)- The 
sequel of this section is devoted to the proof of Theorem [2l The existence 
of solution of non-degenerated model (^^) is splitted in three steps. The 
first one based on approached solutions solving a time discrete system with 
non-degenerate mobilities. For that, let be T > 0, M G N* a number of time 
step, 5t = T/M the time step and let initialize a sequence parametrized by 
St with the initial condition p^. Then, if we consider {p^''^,p^''^) € (L^(f2))^ 
with p[^(p^''^)m(s"''^) > and s"'*^ > at time t„ = n6t, we are searching for 
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a solution {p^~^^'''-,p^~^^''') of the following system 



St 



6t 



- di^{c,{Xrr^'^^D^Vp-^''') = r^^\ 

ty I n+l,e\ n,e 



-r/div(^V(p^+^''- 



r 



n+l 



(4.5) 



where = M^^ + e, with e > 0, with the boundary conditions fl 4.ip . The 
regular izat ion of the mobilities lead to the loss of the positivity on the liquid 
saturation. So, the functions M„ and Z are extended on M by continuity 
outside [0, 1]. 

This technique of semi-discretization method in time has been used by 
Alt and Luckhaus [l| for degenerate parabolic system and has been employed 



m 



151 . llTl . |5[ for a porous medium. A Leray-Schauder's fixed point theorem 
t\ allows to define a solution [p^'^^"" ,p^^^''') for the system ^^ ,5^. 

The second step is devoted to pass to the limit when e goes to zero and to 
prove the positivity of the liquid pressure. A uniform estimate (with respect 
to e) based on the scalar product of fl 4.4p with gg{pg) := /q^" -^^^dz and 

(1 4.5p with ip = Cipi — Pg ensures by using f l 2.2p and f l 3.3p that 

(p^)e, {Pa)t is uniformly bounded in i7p^(n), 

{B{s'i))e is uniformly bounded in H^{^1), 

(Vpc(s;))e is uniformly bounded in L^(f2). 

Up to a subsequence, the sequences (s^)^, (p^)e, {Pa)t, verify the following 
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convergences 



— >p and — Pa weakly in H^^{Q) and a.e. in fl, 

B{s'i) — >-B{si) weakly in H^{Q) and a.e. in Q, 

Z(s'i) — > Z(si) strongly in and a.e. in il. 

Then, pass to the limit as e goes to zero in formulations (1 4.4p - fl 4.5p to obtain 
(Pg^Pi) ^ -f^ri(^) X -^ri(^) solution of 

Jn Jn 

+ [ Kp1ipg)Mi{si)i\/pi-piipi)g)-Vipdx 
Jn 

+ Cl / Kp^{Pg)Mg{Sg) {V P g " P g{p g) g) " V<^dX 

Jn 

+ (Ci-l)7/ / p^{pg)V{pg-pi)-V^dx= I Tg^d: 
Jn Jn 

f ^(^^) ~< ^dx+ / KMi{si){Vpi-piipi)g)-VMx 
Jn ot 

-V V{pg - pi) -Vipdx = / ^^pdx, 
Jn Jn Pi 

pour tout e {H}. 



(4.6) 



(4.7) 



To prove, that the liquid saturation is positive, we consider ijj = (s/)~ in 
f l 4.7p and according to the extension of the mobility of each phase {Mi{si)sl' = 
0) and the fact that (Z(s/)s;~ = 0), we deduce that si > 0. 

The third step is devoted to pass to the limit as 6t goes to zero to prove 
the existence of a solution of the problem (*^3^). For this, we will show some 
uniform estimates with respect to 6t to obtain uniformly bounded on some 
quantities. 

The next lemma gives us some uniform estimates with respect to 6t. 



Lemma 1. (Uniform estimates with respect to St) The solution of (14.4 
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fl4.5p satisfies 

in 



(4.8) 

where C does not depend on 6t. The functions Ti and Vc are defined by 

rsi rPg I 

^giPg) ■= pKP9)9g{Pg) - Pg, 'PM) := / Pciz)dz et gg{pg) = / ^^d^. 

^0 Jo Pi {^) 

Proof. Let forget the exponent n + 1 in this proof and let denote with the 
exponent * the physical quantities at time t„. So, since gg is concave {g^ (p) < 
0), we have 

ggiPg) < 9giP*g) + 9'g{P*g)iPg " P*g) , 

and from the definition of Hg, one gets 

{Pi{Pg)m{si) - p'l{j)*g)m{s*i)) gg{pg) + {si - Si) {Cipi - pg) 
> 'Hg{pg)m{si) - 'Hg{p*g)m{si) - Ci{si - s^)Pc{si). 

Using the concavity of Vc we have the inequality : (s; — s^)pc{si) < Vc{si) — 
Vc{s*i), and the above inequality (14.91) . we obtain the following inequality 

(pf (p,)m(sO - p1{p*g)m{sl)) ggipg) + (s, - s^) {Cipi - Pg) 
>ng{pg)m{si)-ng{pl)m{s\)-C{Pc{si)+CrVM)- 

Now, to obtain the inequality (14.81) . we just have to multiply (I 4.6p by 
9gi.P^^^) (I 4.7p by {Cip^^^ — p^"*"^), sum this two equations and use the 
inequality (1X10]) . □ 

The limit as 5t goes to zero is similar to the limit as e goes to zero with 
additional difficulties on time derivative terms which are overcome in the 



same manner as m 



151,1171,151. 
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5. Existence of solutions of the degenerate system 

We have shown in the previous section HJ the existence of a solution 
iPgyPi) of the problem The aim of this section is to pass to the limit as 
T) goes to the zero to prove the main result of this paper. 

The first point to do this is to obtain uniform energy estimates with re- 
spect to T]. The second point is devoted to gat uniform estimates o, space 
and time translates which provide compactness results on solution by virtue 
of Kolmogorov's theorem. Next, we will be able to pass to the limit as r] goes 
to zero. 

Now, we state the following two lemmas in order to establish uniform 
estimates with respect to r]. 

Lemma 2. The sequences (s^)r, and {p^)r^ satisfy 

almost everywhere in Qt, (5-1) 
is uniformly bounded in L'^(0,T; H^^{il)), (5.2) 
is uniformly bounded in L^(Qt), (5-3) 

is uniformly bounded in L^{Qt)^ a = l^g (5.4) 

is uniformly bounded in L'^{0,T; H^{Q)), (5-5) 

is uniformly bounded in L^(0,T; (ifp^(f2)) ), (5.6) 

is uniformly bounded in L^(0,T; (ifp^(f2)) ). (5.7) 

Proof. The positivity of the saturation fl 5.ip is conserved through the limit 
process. For the next four estimates, we just have to multiply fl 4.2p by 
dgipl) = /o'* p'l(z)Az (I 4.31) by Cipl — and adding them. We follow the 
same calculation as in section [3] to provide the energy estimates fl 5.2p -( r53|l . 

For all if, ^Ij e L^{0, T; H^^{n)) and by using the formulation frT2D - frT3l) 
with the relation fl 2.2p between the pressure of each phase and the global 
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pressure, one gets 



JQt 



+ 



Qt 

+ 



/ XrA'(pf(P^)')VP5-Vv^dxdt + / r.^dxdt 



and 



JQt 

K (Mj(s;')Vp'' + WB{s1)) ■ V^dxdt 



— T/^dxdt 

Qt 

where the bracket (■, ■) represents the duahty product between L^(0, T; {H^^{il))') 
and L^{0,T;H^^{n)). 

From the estimations (15 ■2p and fl5.5p . we deduce 



[<^d,{p^{p^^)m{s'/)),^)\<CM L2(0,T;Hi (Q)), 



and 



which estabhshes f l 5.6p -( r5.7p and proves Lemma O 



□ 



In the next lemma, we derive estimates on differences of space and time 
translates of the function = p'l{p^)m{s^) which imply that the sequence 
(pf (ppm(s^)),, is relatively compact in L^{Qt). 

Lemma 3. {Space and time translate of U). Under the assumptions (HI) — 
{H8), the following inequalities hold : 



\U''{t,x + y) - U''it,x)\dxdt < 



\U'^{t + T,x) -U'^{t,x)\dxdt < u{t) 



(5.8) 
(5.9) 



nx{o,T-T) 
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for all y e and for all r G (0,T); with il' = {x e 9,, x + y e 9,} 
and the function u and Cj are continuous, independent of t] and satisfying 
lim\y\^oUj{\y\) = and \imr^oU}{r) = 0. 

Proof. For the space translates, we observe that 



(o,T)xn' 
{o,T)xn' 



IW^it.x + y) -U'^{t,x)\dxdt 



< 



{o,T)xn' 



P\\p1)m{s^i) )it.x + y)- ip\\pr)m{s'l) ] (t, x) 



m(.n(t, X + y) {p^ip^git, x + y))- pf{p^g{t, x))) 



'{o,T)xn' 

<Si+S2, 



p^{p^„){t,x){m{s'l){t,x + y)-m{s'l){t,x)) 



dxdt 
dxdt 
dxdt 



where £i and £2 defined as follows 



£1 = pM 



(o,r)xn' 



s'^{t,x + y) - 



dxdt. 



(5.10) 



£2 



pKp1{t,x + y))-p\\pm,x)) 



dxdt. 



(5.11) 



'(0,T)xQ' 

To handle with the space translates on saturation, we use the fact that 
B^^ is an Holder function, applying the Cauchy-Schwarz inequality and from 
fl 5.5p . we deduce 



£i<C 



< C 



< c 



{o,T)xn' 

T 



B{snt,x + y))~B{s^{t,x)) 



dxdt 



\/B{s^{t, X + ry)).ydr jdxdt 
Jn' ^Jo ' ■ 

I 1 2 \ 2 

\VB{s'^{t,x + ry))\ drj \y\dxdt 
Jn' "-Jo ' 



<C\y\ 



(5.12) 



To treat the space translates of £2, we use the relationship between the gas 
pressure and the global pressure, namely : Pg = p — p defined in fl 2.2p . then. 
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from the estimation on the global pressure fl 5.21) and the estimate f l 5.12p we 
have 

S2<C{\y\ + \y\'). 

Define = $f/^. From assumption {HI) on the porosity, we deduce thee 
space translates on V^. The proof of the time translates of can be found 
in 111 for more details. □ 



From the previous two lemmas, we deduce the following convergences. 

Lemma 4. (Strong and weak convergences). Up to a subsequence the se- 
quence {s^)rf, {p^ := pI and {pX)r] verify the following convergence 

p^ — >p weakly m L'^{0,T;H^^{n)), (5.13) 

Bis'^) — y B{si) weakly m L\0, T; H\n)), (5.14) 

p'^g — >pg weakly m L'^{0,T;H^^{n)), (5.15) 

— > Si strongly in L'^{Qt), a.e. in Qt, (5.16) 

Si > almost everywhere in Qt, (5-17) 

— > Pa almost everywhere in Qt, (5.18) 

^dtiU'^) ^ ^dt{p'lipg)m{si)) weakly mL\Q,T;Hl^{p)), (5.19) 

^dts^i — y^dtsi weakly m L^{0,T;H^^{n)), (5.20) 

where = p'l{p^)m{s^) . 

Proof. The weak convergences f l 5.13p - fl 5.15P follows from the uniform esti- 
mates f l 5.2p and fl 5.5p of lemma [2J 

By the Riesz-Frechet-Kolmogorov compactness criterion, the relative com- 
pactness of in L^{Qt) is a consequence of Lemma [3] and then ensures the 
following strong convergences 

(f)p'l{p^)m{s^) — > I strongly in L^{Qt) and a.e. in Qt , 

and consequently 

(p^)m(sf) — yU = 1/^ strongly in L^{Qt) and a.e. in Qt , (5.21) 

In order to prove the convergence fl 5.16p . we reproduce the previous Lemma 
|3]for V = (^s^i and as an application of the Riesz-Frechet-Kolmogorov com- 
pactness criterion we establish fl 5.16p . And from f l 5.16p . we deduce fl 5.17p . 
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The convergence fl 5.211) combined with fl 5.16p . to prove that 

p1 — > Pg a.e. in Qt , 

then the convergence flS.lSp for a = is then estabhshed. And again as 
consequence of f l 5.16P with the capillary pressure law, we deduce f lS.lSp for 
a = I. At last, the weak convergence fl 5.19P and f l 5.20p is a consequence of 
the estimate fl5.6p and (15. 7p . 

□ 

In order to achieve the proof of Theorem [H it remains to pass to the 
limit as 77 goes to zero in the formulations fl 4.20 - 4.3p . for all smooth test 
functions and ip in C^([0,T]; H^^{Q)) such that (p{T) = 'ip{T) = 0, 

Jqt 

Kpf (p^MK^n - PiiPi)s) ■ V^dxdt 

JQt 

+ (Ci - l)r/ / pf (pP V(p^ - p1) ■ V^Mt 

JQt 

rgLpdxdt+ / $pf(p°)m(s°)v9(0,x)dxdt, 
Jqt 

(5.22) 

^sldt^dx&t+ I ¥.Mi{sl){Vpl - pi{pi)g) -Vijdxdt 



-V VijPg-p^i) ■Viljdxdt= I ^ipdxdt+ I $s°^(0, x)dxdt, 
JQt JQt Pi JQt 

(5.23) 

The first term in fl 5.22p and f l 5.23p converge due to the strong conver- 
gence of Pi{pl)m{s^) to p^{pg)m{si) in L^{Qt) and the strong convergence 
of to si in L'^{Qt)- 
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The third and fourth term in fl 5.221) can be written as, 




+ / Kp^ipl)VB{sl) ■ Vifdxdt. 



(5.24) 



The two terms on the right hand side of the equation (15.241) converge 
arguing in two steps. Firstly, the Lebsgue theorem and the convergences 



(Tsnej) and (Tsnsj) . estabhsh 

Paipl)Maisl)Vip — > PaiPa)M^{sa)\/ip strougly in (L'^iQT))'', 



Secondly, the weak convergence on global pressure (I 5.13P and the weak con- 
vergence (I 5.14p combined to the above strong convergences allow the con- 
vergence for the terms of the right hand side of (I 5.24p . In the same way for 
the second term of the equation (I 5.23p . 
The fifth term of equations (I5.22p can be written as 



V p'l{p1)V{p'l-p1)V^dxdt = ^ pf(p;')(v^Vpe(5l'))V^da;dt, 



the Cauchy-Schwarz inequality and the uniform estimate (15.31) ensure the 
convergence of this term to zero as r] goes to zero. In the same way for 
the third term of equation (I 5.23p . The other terms converge using (I 5.16p - 
(I 5.18P and the Lebesgue dominated convergence theorem. 
The weak formulations (I 2.7p and (I 2.8p are then established. And The main 
theorem [T] is then established. 
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